Alberta Thy-01-99 
January, 1999 



Long-range two-body final-state interactions and 
direct CP asymmetry in B + —> tt + K° decay 

A. N. Kamal 

Theoretical Physics Institute and Department of Physics, 
University of Alberta, Edmonton, Alberta T6G 2J1, Canada. 



Abstract 

We present a calculation of the direct CP asymmetry, A^jp, for the process 
B + — » tt + K° including the effects of long-range two-body inelastic final-state 
interactions (FSI). We admit four two-body channels in our calculation: B + — > 
(ir + K°), (nK + ), (D^D°) and the charge-exchange channel tt°K + . The strong 
scattering is described in terms of Pomeron and Regge exchanges. We find 
that the direct CP asymmetry is enhanced by a factor of ~ 3 as a result of FSI, 
but remains well short of the claims of (10—20)% in recent literature. A critical 
assessment of papers claiming large CP asymmetries is also presented. Many- 
body decay channels could play an important role in direct CP asymmetry; 
however, their contribution is hard to estimate. 

PACS categories: 11.30.Er,13.25.Hw. 
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I. INTRODUCTION 



Considerable attention [|], |2], |3|] has been paid to B — > nK decays for their poten- 
tial to determine the CKM angle 7. The effect of two-body final-state interactions 
(FSI) has also been studied variously by Refs. [fj, |], ||, §, In calculations |J 
without the inclusion of FSI effects, one finds a direct CP asymmetry, A^p, in 
B + — > 7T + fC° of the order of 0.5%. However, Refs. and calculating the ef- 
fect of an intermediate pK* state converting to tt + K° through an elementary pion 
exchange, estimate the direct CP asymmetry arising out of this triangle graph to 
be ~ 10%. This would imply that such a large CP asymmetry could be accomo- 
dated within the Standard Model. Further, it has been argued in Refs. J7[ that in 
B + — > 7T + K° decay a CP asymmetry of the order of (10 — 20)% is possible once 
long-distance (soft) FSI are included. Hence it is argued that a 20% CP asymmetry 
could not be taken unambiguously to signal New Physics. Ref. || has also dis- 
cussed the role of FSI in the determination of the angle 7, and concluded that CP 
asymmetry as large as 15% can be expected in B + — ► tt + K° decay. 

The goal of the present paper is rather limited: it is to investigate the effect 
of long-range two-body FSI on Aq P in the process B + — > tt + K° which in the ab- 
sence of FSI is a pure penguin- driven process. We admit four two-body channels in 
the discussion of inelastic final-state interactions: (n + K°), (r]K + ), (D+D°) and the 
charge-exchange (CEX) channel tt°K + . The decay amplitude is decomposed into 
the two isospin states, I = 1/2 and 3/2, and each corrected for FSI. Our method 
of including FSI is different from that of [[/J. In contrast to the claims in Ref. 0, 
we find that A^ P could be raised by a factor of ~ 3, but remains well short of (10 
- 20)% range. We spell out the role played by the elastic and inelastic two-body 
channels on Aq P , and take a critical look at Refs. |5|, [F| thereby tracing also the 
reasons for the variance between our results and theirs. Following Ref. we also 
neglect the electromagnetic penguin processes. 

The layout of the paper is as follows: We define the Hamiltonian and the model 
parameters in Section II. Weak decay amplitudes for the four two-body channels, 
B+ _> (n + K°),(r]K + ), (D+D°) and (tt°K+), are calculated in Section III. The 
strong S'-matrix elements, Su, S12, 5*13 and Su, are calculated in Section IV, using 
Pomeron and Regge exchanges. CP asymmetry calculations are presented in Section 
V. The concluding section, Section VI, contains a discussion of the results, and a 
critique of previous calculations claiming large asymmetries. 



II. DEFINITIONS AND PARAMETERS 
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In our analysis we will consider only the the QCD-modified effective Hamiltonian 



for b — > s transitions with the top quark integrated out [1C, 11, 12 



Hejf = % £ \y qh y; s \c x o\ + c 2 oi + £aoA (i) 

V ^ q=u,c V j=3 J 

Here Cj are the Wilson coefficients defined at scale /x ~ and the current 
bilinears O, are defined as follows: 



0? = (sq)v-A(qb)v-A, 

3 = (sb) v ^ A J2 q '(q'q')v-A, 



5 = {sb) V - A Y.q< {q'q') 



V+A, 



oi 
o 4 



{s a q/3)v-A 

{Sah)v-Aj2Wo 



V-A, 
V-A, 



(s a bp) V - A ^2(^ a )v+A- (2) 
g' 



In Eq. (2), 0\ and 2 are the tree operators, 3 , ....,0 6 are generated by QCD 
penguin processes, (V ±A) represents 7 M (1 ±75), and a and (5 are color indices. J2 q ' 
is a sum over the active flavors u,d,s and c quarks. 

In the next-to-leading-logarithmic calculation one works with effective Wilson 
coefficients C\ , rather than the coefficients that appear in (1). The derivation of 
these effective coefficients is well known [Rll, O, [12|. We simply quote their values 



u 1 — Oi, u 2 — O2, o 3 — O3 — r B /i\ c , o 4 — O4 + -r s 



withiTl 



(3) 



and 



where 



Ci = 1.1502, C 2 = -0.3125, C 3 = 0.0174, 
C A = -0.0373, C 5 = 0.0104, C 6 = -0.0459, 



-f^Ci(//) [- + - In -f - G(m„ //, g 2 )] , 

o7T 3 O jJL 



G(m„, /1, g 2 ) = — 4 / c/xx(l — x) In [1 — x(l — x) — -]. 

Jo mi 



(4) 



(5) 



(6) 
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Note that only the current operator Of, with q=u,c, is used to generate the penguin 
contributions; q 2 , on the other hand, is the invariant momentum carried by the gluon 
in the penguin diagram, and m q the mass of the quark q in the penguin loop. For 
q 2 > 4m 2 , G(m q , fi,q 2 ) becomes complex giving rise to strong perturbative (short 
distance) phases through P s . We use N c = 3 throughout. 

The quark masses and the meson decay constants, we employ are: 

m u = m d = lOMeV,m s = 175MeV,m c = lMGeV,m b = 5.0GW, 

f w = UlMeV, Ik = 160MeV, f v = UOMeV, f Ds = 300MeV. (7) 

As for the CRM angles, Ref. [0 quotes \ V ub \ / \ V cb \= 0.08 ± 0.02, and Ref. 
[0 quotes Arg (V* b ) = (62 ± 12)°. Refs. || and pf cite a significantly different 
range for Arg (V* b ). For our calculations, we take, 



V, 



ub 



cl, 



0.1, and Arg(V:j = (50 o ,75 o ,90 o ). (8) 



Since long-range rescattering preserves isospin, we decompose B + — > tt + K° am- 
plitude into 1=1/2 and 1 = 3/2 components and then apply Pomeron and Regge 
exchanges to these amplitudes to compute the FSI-corrected isospin amplitudes. 
Finally the FSI-corrected isospin amplitudes are recombined to generate the total 
FSI-corrected decay amplitude which is then used to calculate A d Q T P . 



III. DECAY AMPLITUDES 

A. Channel B+ -> n+K° 

The B — > tiK decay amplitudes in terms of the isospin amplitudes are given by 
(see, for example, Ref. 0]) 

A +0 = A{B+ ^n+K ) = A 3/2 + A ( Jl 
A)+ = V2A(B + -> n°K + ) = 2A 3/2 - A^ 2 
A _ + = A(B° ^n-K+)=A 3/2 + A{-l 

A 00 = V2A(B°^n K°) = 2A 3/2 -A[-/l (9) 

where 

Ayl = A 1/2 + B l/2 , A { -\ = A l/2 - B 1/2 , (10) 
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and 

A 1/2 = ± v /2/3(l/2,±l/2|iJ A/=1 |l/2,±l/2) 

B 1/2 = ^73(l/2,±l/2|// A7=0 |l/2,±l/2) 

B 3/2 = V / l73<3/2, ±l/2|i? Ajr=1 |l/2, ±1/2). (11) 

In Eq. (11) the state vectors are eigenstates of / and I 3 . if a/=o and i?Ai=i are 
the isospin-conserving and isospin-changing parts of the Hamiltonian, respectively. 
Note that the penguin operators, O3 to 0$, involving ^2 q (qq) are of necessity isospin 
conserving. The operators Oi and 2 generate both AI = and AI = 1 transitions. 
Note also that there are three independent isospin amplitudes but four decay am- 
plitudes. The constraint is provided by requiring that A 3 / 2 calculated from A +0 and 
Aq+ be the same as that from A- + and Aqq. This serves to fix the relative signs of 
the decay amplitudes, as we shall see in the following. 

In the factorization approximation the four decay amplitudes of Eq. (9) are given 
by ( an overall factor of Gf/V% has been suppressed) 

A +0 = -fK{ml-ml)F^{ml){aa { : ) + R^a^) 

+C(4 c) + *f ) 4 c) )}, 

A 0+ = fK(ml-ml)Ft(ml){C u (a 1 + ai u) +R { d K) a^) 

+C (4 C) + Rf } 4 C) )} + Um% - m 2 K )F BK (ml)£a 2 , 

A- + = f K (m% - ml)Ft{ml) (a, + a?> + R[ K >a<«>) 

+C(a ( ; ) +Rl K) a^)}, 
Aoo = -fK{m 2 B -ml)Ft{m\){^ (a™ + R d K) a^) 

+C (4 C) + Rf } 4 C) )} + U(m% - m 2 K )F BK (ml)£a 2 , (12) 

where 

«2i — <^2i + ^~°2i-l> 

o«_i = C^U + ^-C e J f , (N c = 3,i = 1, ...,6) 

ft(K) _ 2m i~ ^ _ v b V s * (13) 

q (m b - m q )(m s + m q y q qb qs ' 
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and the superscript on af^ and represents the quark flavor in the penguin loop. 

The isospin amplitudes derived from the above equations are (note that one has 
to set = Rf ] for isospin symmetry) 



,<*)<>) 



M/2 = y:{(ml-ml)f K Ft(ml)a 1 + (ml-ml)UF BK (ml)a 2 } 

4^ = -le u {(rn%-ml)f K Ft(m 2 K )(a l + 3a { t ) + 3R^ 
+ (m 2 B -m 2 K )f 7r F BK (ml)a 2 } 
-€« - ml)f K Ft(m 2 K ) (4 C) + Rf^f 



472 = 56 {(m|, - ml)f K Ft(m 2 K ) (201 + 3af + 3^ 
-(m| - mD/^K)^} 

+€(™l - O/xi^ «) (4 C) + J#°a< c >) • (14) 



Notice that A3/2 gets contribution from the operators 0\ and 2 only. For the 
discussion of B + — > ty + K° decay only A 3 / 2 and Ayl are relevant. With the chosen 
parameters of Eq. (7), and the form factors of fl7fl , they are (aside from a factor 

G F \. 

a 3/ 2 = c4 3/2) + C4 3/2) 

= ^(0.5529) GeV 3 
a(+) = t* T (i/2) , (1/2) 

^1/2 — Sn^u T Sc^c 

= £(-0.4637 + 0.04127i)+ £(0.1034 + 0.03595i) GeV 3 . (15) 

Note that x^ 3/2) = 0. The two isospin amplitudes are unitarized through isospin 
conserving FSI. Henceforth, the superscript (+) in A^l will be dropped wherever 
it is possible to do so without ambiguity, since it is the only relevant isospin 1/2 
amplitude. 



B. Inelastic Channels B + -> r]K + , B + -> D+D° and n°K + 

We specifically consider the effect of three channels, B + — > rjK + , -B + — > DfD° 
and the charge-exchange channel S + — > 7r°if + , on the elastic channel S + — > vr+i^ . 
The first two inelastic channels will effect only Ai/ 2 , while the third will effect 
both A1/2 and A3/2. In the following we enumerate the three decay amplitudes 
B + -»■ D+D , 5+ -» and B + -> vr ^. 
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r]K + Amplitude 



Due to i] — rj' mixing, the expression for the decay amplitude for B + — > rjK + is 
fairly involved. In writing the following m 2 B — m 2 K « m 2 B — m 2 « has been used 
( a factor of ^§ is suppressed) , 



04 [^^(O^(^) (oi + of + r { u k) c 

+Fo BK H)f v [F(9)a 2 + ^ (af° - a™) 

+G(9) (aP + RMaP)}] 

H> B 2 [f K F Q B \m 2 K )F{e) (ai c) + R^af 

,rM/ 2 U f sin6 /(«) fu)\ 

+G(0) (a?> + i#>a<">)} 



(u) 



(16) 



where 



2m^ 



(m 6 - m s ){2m s y 



G(6) = - ( 4=cos# + * sin # 



6 = 7] — t]' mixing angle. 



(17) 



The largest contribution to the decay amplitude arises from the tree amplitudes 
proportional to a\ and a 2 which constitute the coefficient of £*. Numerically, with 
9 = —20°, and form factors from IfC 



V K + ) = CVu + CVc 

= £*(0.8901 - 0.00318i) 

+£(-0.640 - 0.277i)10" 2 GeV 3 . 



is) 



B + -> D+D° Amplitude 
In the factorization approximation one obtains ( a factor of is suppressed) 

A(B+^DfD°) = (m|-m^)/ D .F^(m!,j{c(afJ + iJ< D ->4 tt) 

+c(«i + 4 c) + ^ s) 4 c) )} 



(19) 
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where 

R? a) = 7 ^ \ (20) 

(m b - m c )(m s + m c ) 

Numerically, 

a{b+ - DjD°) = ex+c^ 

= £*(-0.4758 - 0.2067i) 

+£*(4.5055 - 0.1800i)GW 3 (21) 

We have used fn s = 0.3 GeV and F BD (mD 2 ) =0.66 from Heavy Quark Effective 
Theory 



B+ -> 7r°^+ Amplitude 

From Eq. (9) 

where A 3 / 2 and A^ have been defined in Eq. (15). 



A(B+ -> vr ^) = - -U« (22) 



IV. INELASTIC FINAL-STATE INTERACTIONS 



A. Formalism 

For the purposes of the following discussion, label channels ir + K°, DfD°, rjK + 
and tt°K + as channels 1, 2, 3 and 4, respectively. Let us also denote by the matrix 
&ij (hJ = 1) 2, 3, 4) the elements of the 5-wave scattering matrix in isospin / state. 
Then the effect of rescattering on the weak decay amplitudes Aj is to generate the 
unitarized amplitudes via |0 



A u m = ( i + *>) (23) 
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where A- 7 "* are the weak decay amplitudes calculated so far in this paper according 
to the isospins. For the two isospin amplitudes, the above equation reads explic- 
itly (note that the isospin label is now a subscript in the decay amplitudes but a 



S 



superscript in the S matrix elements), 

A u 1/2 = l(l + s[^-^Sir^A 1/2 (B^^ 

/ 0(1/2) \ 




^13 



(1/2) 



^vK + ), (24) 



where we have used Eq. (22) in writing the contribution of channel 4. The am- 
plitudes A 1/2 (B + — >• n + K°) and y4 3 / 2 (i? + -> 7r+^°) are given in Eq. (15), and 
A(B + -> ?7J^ + ) and A(5+ -> £>+l) ) in Eqns. (18) and (21), respectively. The 
effect of the charge-exchange channel tt°K + has been simply to affect the following 
replacements: 

o(l/2) _^ o(l/2) _ 1 Q (l/2) 

°11 ~~ * °11 ' 

S§!*> - sl 3 / 2) + V2S{T } ■ (25) 

Our next task is to calculate the scattering-matrix elements. For this we use the 
Regge-exchange model including the Pomeron and the highest lying trajectories for 
each reaction. The S matrix for the scattering from channel i to channel j in L = 



state is given by |L9 



Sij = 5ij + _ / T(s, t)dt, (26) 

87Ty AjAj Jtmin 

where t m i n and t max are the minimum and maximum values of the invariant mo- 
mentum transfer t, and Aj and Xj are the usual triangular functions: X(x,y,z) = 

(x 2 + y 2 + z 2 — 2xy — 2xz — 2yz) 1 ^ 2 for channels i and j, respectively. For u channel 
exchanges, replace t by u. In the following we consider the evaluation of S^. 



B. Evaluation of S 



n 



For the elastic scattering, n + K° — > tt + K°, we consider only the Pomeron, the p and 
/ trajectories. K* exchange in the u channel is not permitted in tt + K° — ► tt + K° 
scattering as the exchanged object is an exotic (qqqq) state. A (qq) exchange is, 
however, allowed in n + K° -> tt + K° scattering;. 
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Pomeron Exchange 

The elastic scattering 7r + K° — > 7r + fr° via Pomeron exchange is parametrized by [[20 



Tp(s,t)=P(t)(-) J™ p (i)/2 , (27) 

with |2T| 

a P (t) = 1.08 + 0.25t, (28) 



and 21 



0(t) = /3(0)e 2 - 84 , (29) 

where the invariant momentum transfer t is expressed in units of GeV 2 and its 
coefficient in units of GeV~ 2 . The energy scale, Sq is taken to be 1 GeV 2 for 
7i + K° — > n + K° scattering. The coupling constant /3(0) is fixed by the additive 
quark model, yielding /3(0) = 2/3 un + 2/5 us , where the flavor subindices refer to the 



quark flavors between which the Pomeron is exchanged. Using [3 US = \(3 UU [p2 



suggested by high energy irp and Kp elastic scattering data, and (3 UU ~ 6.5, we 
obtain /3(0) ~ 22. With these parameters we calculate, 

[ T P (s,t)dt = -0.313. (30) 
87tAh J 

p Exchange 

The p exchange amplitude for n Q K — > tt^K scattering is parametrized as 

TW\s,t) = ±[T P ,T a ]T p (s,t), (31) 

where a and (3 are the charge indices of the initial and final state pions, respectively, 
and 

T P (s,t)=P p (0) } Lr(-) , (32) 

r{a>p{t))smira>p(t) \s J 

where s is again chosen to be 1 GeV 2 , and the p trajectory is taken to be, 

a p (t) = a p (0) + a' p (t), (33) 
with a p (0) = 0.5 and a' = 0.848 GeV' 2 . 
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The isospin projection operators in I — 1/2 and 1 = 3/2 states are, 



Qa/> 2) = g 7 /^ Q% 2) = <W ~ ^r p r a . (34) 

Thus, 

\[r,,r a }=2Q^-Q% 2 \ (35) 

and p exchange contributes to the two isospin states / = 1/2 and I = 3/2 in 
the ratio 2 : — 1. The coupling constant (3 p (0) is determined by taking the limit 
(s — > oo,t — > m^) in Eqns. (31) and (32) for the reaction tt+T^ — > 7r + fT and 
equating it to the perturbative graph contribution with the p meson in the t channel. 
This leads to 

0,(0) = -^9 2 vpp = 25.1, (36) 

where g VPP /47r = 3.0, determined from the p width, has been used. With these 
parameters and approximating r(a p (t))sin7ra p (t) r(a p (0))sin7ra:p(0) = y/n, one 
gets, 

— ^— / T p (s, i)dt = -0.0202 + 0.01911 (37) 
87rAi J 

f Exchange 



The /-exchange amplitude for n a K npK scattering is given by 

TW\s,t) = 6 a pT f (s,t), (38) 

where, 

^''' Kw^) u) • (39) 

The scale factor s is chosen to be 1 GeV 2 . Since 5 a p = + <3a/ 2 \ f exchange 

contributes in equal measure to both isospin states. We consider the / trajectory 
to be parallel to the p trajectory, a'f = a' p = 0.848 GeV~ 2 and determine ct/(0) 
from the central value of / meson mass, 1.275 GeV. This yields «/(0) = 0.6215. 
We also considered a / trajectory degenerate with the p trajectory, but the effect 
of such a trajectory on the CP asymmetry was not significantly different from that 
of non-dgenerate trajectory. Consequently, we present here only the discussion per- 
taining to a non-degenerate / trajectory, albeit parallel to the p trajectory. Further, 
we approximate T (a f(t)) smn a f(t) ~ r (a/(0)) sin7ro;/(0) = 1.3381. The constant 
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/3/(0) is determined by taking the limit (s — > oo,£ — > m/ 2 ) in Eq. (38) (and (39)) 
and equating it to the perturbative expression for the tt + K° — > 7f + K° scattering 
through / exchange. The coupling of / to pion pair and KK is determined from 



the branching ratios B(f — > nn) and B(f — > fCfT) given in [13 . This results in 
/3/(0) = —77.21. With these parameters, we get 

— ^— [ T f (s, t)dt = -0.1577 - 0.29521 (40) 
87rAi J 

Using crossing matrices, the element Su in I — 1/2 and 1 = 3/2 states is given 

by 

S[\ /2) = 1 + [{T P (s,t)+ T f (s,t)+2T p (s, t)}dt 

S{T ] = 1 + [{T P (s,t)+T f (8,t)-T p (s,t)}dt, (41) 

where the subindices refer to the exchanged trajectory. 

Putting the contributions of the Pomeron, the p and the / from Eqns. (30), (37) 
and (40) in Eq. (41), we obtain 

S[\ /2) = 0.4889 - 0.25702, S{l /2) = 0.5495 - 0.31431 (42) 



C. Evaluation of S 
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The inelastic process tt + K° — > r]K + , a pure 1 = 1/2 reaction, proceeds via the 
exchange of K* in the u channel. The Regge amplitude for this process is given by 

T K *(s,u)=p K *(0) K . Vr - • ( 43 ) 

1 {aK*{u))sm7raK*{u) \s J 

The coupling constant /3^*(0) is determined by taking the limit (s — > oo, u — » m* K 2 ) 
in Eq. (43) and equating it to the perturbative amplitude for 7r + i^° — ► r]i^ + process 
through K* exchange. Relating K* —K — rj and K* —K — it couplings to the coupling 
gvpp through SU(3) symmetry, we get 

Pk*{0) = — =gvppna K *'socos9, (44) 

where 9, the rj — rf mixing angle, is taken to be —20°. The scale parameter sq 
is chosen to be 1 GeV 2 and the K* trajectory is chosen to be parallel to the p 
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trajectory (a' K * = 0.848 GeV~ 2 ). This requires a^»(0) = 0.325 for the K* mass to 
be 895 MeV. With these parameters we obtain /3k*{0) = 57.78. In projecting the 
L = partial wave according to Eq. (26), we use the approximation, 

r(a K *(t))simra K *(t) w T (a K *(0)) sin7ra x *(0) = 2.345. (45) 

We finally obtain 

S12 = q — 7=== / T K *(s,u)du 

07TV AlA2 J «min 

= (-0.70448 + 1.10972i)10~ 2 . (46) 

D. Evaluation of Si 3 

The reaction 7r + i^° — > D s + D° is a pure I = 1/2 process. It proceeds through the 
exchange of D* in the t channel. (Here t channel is defined by the reaction tt + D° — > 
K°D S + ; note also that D s * cannot be exchanged in the u channel, tt + D s ~ — > K°D° 
since the exchanged object has to be neutral.) The Regge amplitude is 

T 13 ( s ,t) = (3 D *(o) v -m(-) ■ (47) 

r(a D *(t))smira D *(t) \s J 

Now, the reaction n + K° — > D s + D° has a threshold close to 4 GeV. Thus, we antic- 
ipate the scale factor so to be much larger than 1 GeV 2 , as has been assumed so far. 
Further, if one puts D* (J p = 1 + ) and D2* (J p = 2 + ) on the same linear trajectory, 
one gets a much flatter trajectory, a' D , ~ 0.5 GeV~ 2 , than the p trajectory [|l^, |23[ . 



With this slope one finds a^* (0) = —1.02. In contrast, had we chosen the D* tra- 
jectory to be parallel to the p trajectory, we would have obtained a£>*(0) = —2.426, 
a much lower intercept and, consequently, a much smaller contribution to T\ 3 (s,t). 
In addition, we use the approximation 

r(a D *(t))sm7ra D *(t) rs r(a jD *(0))sin7r« D .(0) = 7r, (48) 

which is due to the cancellation of the pole of the Gamma function at argument 
— 1 with the zero of the sine function at «£)*(0) = —1. The constant (3d* (0) is 
obtained by taking the limit (s — > 00, t — > m 2 D *) in Eq. (47) and equating it to 
the perturbative graph with D* exchange. The coupling (D*Dti) is related to gypp 
through SU(4) symmetry. This results in /?d*(0) = 118.44. With these parameters, 
we obtain 

sf 3 /2) = - ' I*™ T D *(s,t)dt 



-0.22834 + 0.36992i)10~ 2 . (49) 
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E. Evaluation of Su 

Consider now the charge-exchange (CEX) reaction ir a K° — > tt /3 K + . It proceeds 
through the exchange of the p trajectory in the t channel and that of the K* trajec- 
tory in the u channel. Let us write these two contributions in the following form: 

p trajectory: ^fo, T a ]T^ EX (s, t), 



K* trajectory: T a T p T^ EX {s,u), (50) 



where 



l-e" im ' (i) ) / S W*) 



T(a p (t))sm7ra p (t) Vso 



T^j (s,u) = (3 { c K EX (0 -A , Vt(-) • 51 

r(a^»(-u))sm7ro;ft:.(-u) \s / 

The constants (3 { C P EX (0) and /$S(0) are determined as in other cases heretofore, to 
be 

PclxiO) = 2/3,(0) = 50.21, /Ci(0) = -0 P (O) = -25.1. (52) 
Using the crossing matrices for the isospin spin operators, we obtain 

S{\ /2) = -!— (2 f maX T^ EX {s,t)dt - r aX T^; x (s,u)du) 

07TAl \ Jt min JUmin ' 

= -0.08386 + 0.081207i (53) 

S[ 3 J 2) = -1- (- f maX T^ EX (s, t)dt + 2 I"™* TS( S , «)d«) 

= 0.046527 - 0.047815i. (54) 



V. CP ASYMMETRY CALCULATIONS 



A. Definitions and direct CP asymmetry without long-range FSI 
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Let us define the generic form of the decay amplitude for B + — > tt + K° as 

A(B + -> n + K°) = C u Au + £A C , (55) 
then the decay amplitude for the process B~ — > tt~K° is 

A(B- -> n-R°) = i u A u + £ c Ac (56) 
We define the direct CP asymmetry as, 

A * = - k+K°) - T(B- -> vr-iT ) 

CP r(s+ ^ tt+x ) + r(5- ^7T-i? )' 1 J 

Further, if we define the CKM products with their weak phases following Wolfen- 
stein parametrization, 

C u =\iu\e l \ 6H&|, (58) 
and the decay amplitudes A u and A c with their strong phases, 

A u = | A u | e ", v4 c =| A c | e c , (59) 



then 



where we have defined, 



A dir „ -2r£sm(5 u - £ c )sin 7 

^^1 + 2^08(^-4)0087' 1 J 



VubV us Vub j / e -i\ 

(= \T\ = Twi\ =x \v*\- aair= \ATV (6l) 

where A = 0.22 is the usual Wolfenstein parameter representing V us , and in the 
denominator of Eq. (60), we have neglected a term of order | £ | 2 . Note also from 
the expression for Afj r p , Eq. (60) (which corresponds to the definition used in Ref. 
[24|j), that for a given value of the angle 7 the CP asymmetry can be raised either 
by raising r, or by raising sin(5„ — 5 C ), or both. This point is discussed in the next 
section in the context of our calculations. 

In the absence of long-range FSI, the decay amplitude for B + — > tt + K° is given 
by A + o of Eq. (12). One can then calculate the direct CP asymmetry from the 
numerical values for A 3 / 2 and Ai/ 2 given in Eq. (15) with A +0 = A 3 / 2 + Ayi. Recall 
that the superscript (+) in A^l has been dropped. 
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With A = 0.22 and £ = 0.022 (see Eqns. (8) and (61)) and the three choices for 
7 = (50°, 75°, 90°), the CP asymmetries calculated with q 2 = m 2 /2 are (and shown 
in Table 1) is, 

Af P = (-0.29%, -0.37%, -0.39%), for 7 = (50°, 75°, 90°). (62) 

The reason for the asymmetry being an order of magnitude smaller than | £ |~ 0(\ 2 ) 
is that sin(5 M — 5 C ) is of order 10 _1 (see Table 1). Note also that the asymmetries 
in absence of FSI are comparable to those calculated in M. The authors of Ref. || 
have used a different set of mass parameters; the phase 5 C is sensitive to the choice 
of m c , particularly when m 2 c is in the vicinity of m 2 /8. In fact, the CP asymmetry 
is quite sensitive to the choice of q 2 (see Ref. 0). The difference in the sign of the 
asymmetry between us and Ref. is a matter of definition only. 

B. CP asymmetry with long-range inelastic FSI 

In terms of the two isospin amplitudes, A 3 / 2 and Ayl, defined in Eq. (14), the 
decay amplitude for B + — > tt + K° is given by 

A +0 = A(B + - n+K°) = A 3/2 + A$. (63) 

After unitarization, the two relevant isospin amplitudes are given by Eq. (24). 
Substituting from Eq. (24) in Eq. (63) and reorganizing the resulting expression, 
we obtain the decay amplitude after the inclusion of long-range inelasic FSI as, 

A +0 = C U A U U + C C A U C) (64) 

where 




The weak decay parameters (x„ , , y u , z u ) and (x£' 2 \ y c , z c ) are given in 
Eqns. (15), (18) and (21). The S-matrix elements , , S12, Si 3 , S^ 2 * 1 and 
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S lA are given in Eqns. (42), (46), (49), (53) and (54), respectively. The CP 
asymmetry is now given by 

A dir _ -2r'£sm(5' u - ^)sin 7 
CP ~ 1 + 2r'£cos(<^ - ^)cos 7 ' 1 ' 

where 

A%=\AZ\e< A u c =\A u c \e i5 ' % r> = j^j. (67) 
The resulting CP asymmetries calculated with q 2 = m 2 /2 are: 

t4*p = (0.73%, 0.94%, 0.99%), for 7 = (50°, 75°, 90°) (68) 

Comparison with the CP asymmetries in absence of FSI shows that CP asym- 
metries have increased approximately three-fold, and have changed signs. 

We studied another scenario - one where channels 2 and 3 were removed and 
only channels 1 and 4 retained. In such a scenario only the elastic and the CEX 
channels come into play and the effect of the FSI are embodied in S$ and S$- 
The results were, 

Afp = (0.50%, 0.65%, 0.68%), for 7 = (50°, 75°, 90°) (69) 

We note from Eqns. (68) and (69) that adding channels 2 and 3 has made a signif- 
icant difference to the asymmetry. 

In yet another scenario in which only channel 1 was kept (only elastic re-scattering 
allowed), the CP asymmetries were, 

Af p = (0.44%, 0.56%, 0.59%), for 7 = (50°, 75°, 90°). (70) 

The results are summarized in Table 1. 



VI. DISCUSSION 



In the absence of FSI, the decay B + — > ix + K° is a pure penguin process; the decay 
amplitude depends on the coefficients Oj(i = 3, 6). With our choice of the parame- 
ters, and with an effective q 2 = m 2 /2, we have calculated the direct CP asymmetry 
A d Q P in the absence of FSI to be in the range (—0.3 to — 0.4)%, depending on the 
choice of 7. This is an order of magnitude smaller than A 2 . 
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Once the decay amplitude for B + — > tt + K° is decomposed into the two isospin 
states, 1 = 1/2 and I = 3/2, each isospin amplitude involves the tree coefficients a\ 
and 02; it is just that they are cancelled in the sum. As the rescattering S matrices 
and are different, the cancellation of the coefficients ai and 02 does not 

occur after FSI are taken into account. That is, FSI mix and blur the classification 
of processes by their topologies, tree or penguin. 

Considerable attention has been paid || || ]j| to the size of the direct CP asym- 
metry in B + — > 7T + K°. It has been argued in that asymmetries as large as 20% 
could be generated as a result of FSI. Refs. [|, |(| argue that 10% asymmetry could 
be produced in the mechanism they propose. Our calculations suggest otherwise. 



We expand on a critique of these papers, and a related one |25| , in the following. 

In order to discuss the results of J7|, we need to establish the correspondence 
between their notation and ours: 

Before FSI: | r + | (Ref.[7]) = | r£ | (here 



VubVus* 




A u 


1 v cb v cs * 




A c \ 



1 v ub 




A u 


1 v cb 




A c \ 



(71) 

Consequently, since | A u \ / \ A c |~ 1, see Table 1, and | V ub \ / \ V cb \ = 0.08±0.02 
[|l]|, one finds | r + |< A/10. After FSI are taken into account, the correspondence 
reads, 

After FSI: | e | (Ref.[7]) = A \ ^ b | [ , ^" | (here). (72) 



It is argued in that e could be as large as 0.09 as a result of FSI. For this to 
be achieved one would need | Aj \ / \ A c ' |> 4. Our calculations show (see Table 
1) that after incorporating inelastic FSI, this ratio has remained close to unity. 
As a consequence, on this score alone, our calculated CP asymmetry would be a 
factor of four lower. A further lowering (by a factor of five) comes about as follows: 
Table 1 shows that the strong phase factor sin(5 u — 5 C ), achieves a value of ~ 0.2 in 
magnitude. In contrast, this factor is not estimated in Ref. [0, and an optimistic 
value of unity is assumed for this factor. This provides a further decrease in the CP 
asymmetry by a factor of five. Thus, our calculated CP asymmetry is a factor of 
20 lower than that in Ref. J7|. 

The estimates of 0, which appear to be too optimistic, could be criticized on 
several grounds. First, to estimate the effect of FSI, they calculate the discontinuity 
of the decay amplitude across the unitarity cut using a Regge amplitude for the scat- 
tering amplitude. A dispersion relation is then used to generate the FSI-corrected 
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amplitude. Their final result (Eqn. (2.17) of Ref. [0]) amounts to retaining only the 
absorptive part of the integral (pole contribution) and neglecting the dispersive part 
(coming from the principal part integration). Besides, the integration ranges from 
s = (m K + m n ) 2 to infinity; Regge description hardly applies at the lower end of 
the integration range. Secondly, having estimated the effect of 7r°K + state through 
the exchange of the p trajectory on the parameter e (see Eq. (2.25) of 0), the con- 
tributions of other channels are estimated only very crudely (see material following 
Eq. (2.25) of [0]). Our estimates do not support the optimistic expectations that 
inelastic channels add up coherently to generate an e of the size required to produce 
a (10 - 20)% CP asymmetry. 

An estimate of a 10% CP-asymmetry in B + — > tt + K° in Refs. (H,f|) is based 
on the calculation of the triangle graph representing the decay B + — > tt + K° via an 
intermediate pK* state. The transition pK* — > tt + K° is assumed to occur through 
an elementary ir exchange. This is tantamount to saying that one-pion-exchange 
mechanism reasonably describes the reaction pK* — > tt + K° at ^/s w 5 GeV . That 
this is not so is discussed below. 

Consider the scattering of helicity zero p + and K*° into tt + K°. The scattering 
T-matrix for this reaction with the exchange of a 7r° is given by 

T ( s > *) = wf P 2 > g * - kE p cos9)(pE K - kE K .oos0), (73) 

m p m* K (t — m~J 

where t = (p n — p p ) 2 = m 2 p + m% — 2E p E n + 2pkcos8. Here 6 is the scattering angle 
in the center of mass, and k and p are the center of mass momenta in the ir + K° and 
p + K*° channels, respectively. 

For zero-helicity scattering, the rotation matrices d,Q (8) are the same as the 
Legendre polynomials Pj(cos6). Hence the projection of J = partial wave proceeds 
as in Eq. (26). Converting the variables from the invariant momentum transfer t to 
cos9, one gets 

S J=0 (pK*° -> n + K°) = [ +1 T(s,cos9)dcos9 

8n^s J-i 

w -28z. (74) 

Unitarity requires that the off-diagonal (and the diagonal) elements of the S matrix 
be bounded by unity. Thus, the one-pion-exchange mechanism leads to a gross 
violation of unitarity of the S matrix. It is for this reason that a large CP asymmetry 
is generated in |JJ and 0. 

Finally, we would like to comment on a related work of Ref. [p5|| . In this paper, 
the authors evaluate the isospin 1/2 and 3/2 scattering amplitudes in nK —> irK; 
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the mechanism is assumed to be through the exchanges of the Pomeron and Regge 
(p, K*, etc.) trajectories. The partial wave amplitudes are then projected and the 
scattering phases, 5i identified as (see following Eq. (26) of ]2"5| ): 



ImAi(s) 
ReAi(s) 

This identification is correct only if the scattering is elastic, i.e. the S matrix is 
given by e 2tSl . To wit: 

If Si = e 2iSl 

then Ai ee i(Si-l) 

= sm5ie iS \ (76) 

where Ai is the scattering matrix constrained by unitarity of the S matrix. Clearly 
now, Eq. (75) holds. However, if the scattering is inelastic, as we anticipate it to be 
in the S-mass region, then each diagonal element (ii) of the S matrix can be written 

as, 

4 = ^, (77) 
where the elasticity rji is such that < rji < 1. Now, 

1 



1 



2 .(^-l), (78) 

and one no longer obtains Eq. (75). Thus, we believe, the identification of the 
scattering phase shift in |25] is incorrect. However, using the real and the imaginary 
parts of the calculated amplitude A\, one can calculate both rji and 8\. 

So much about the scattering amplitude and its phase, the phase of the weak 
decay amplitude is quite a different matter. Consider a spinless single channel 
(elastic) scattering, and a single weak decay amplitude, A(s). The 5* matrix in 
L = state is S = e 2tS . Then the unitarized weak decay amplitude A u (s) is, 

A u (s) = C-^A)A(s) 



2 

= cos5(s)A(s)e iS{s) . (79) 

Clearly, the phase of the weak decay amplitude is the scattering phase (Watson's 
theorem) if the decay amplitude is real. The claim of |25] that the phase they are 
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calculating is the phase of the weak-decay amplitude, therefore, is correct only if the 
scattering is elastic. The scattering, however, is not elastic as evidenced by | Su | < 1 
in our calculation. 

In general, for inelastic scattering, the unitarized weak decay amplitude is a 
superposition of the form given in Eq. (23), and the phase of AV is, obviously, not 
simply related to the phases of the ^-matrix elements, but to the ^-matrix elements 
weighted by the weak-decay amplitudes A,-. 

A final word of caution: Our calculations have been done within the ambit of a 
few, hopefully important, two-body channels. However, our results contradict the 
conclusions of other calculations with the same limitations. Many-body decay chan- 
nels could, in principle, play an important role in determining the CP asymmetry; 
yet their effect is difficult to estimate. Thus, whether large CP asymmetries could 
be generated within the Standard Model remains an open question. 

An Individual Operating grant from the Natural Sciences and Engineering Re- 
search Council of Canada partially supported this research. 
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Table 1: The ratio r =| A u \ / \ A c \, the strong phase factor sin(<5 u — 5 C ), and CP 
asymmetries in % in different scenarios. 



Scenario 


r 

= l ^« 1 / 1 | 


sm(5 u - S c ) 


CP asymmetry 


7 = 50° 


7 = 75° 


7 = 90° 


No FSI 


0.898 


0.098 


-0.29 


-0.37 


-0.39 


Elastic FSI 


1.071 


-0.125 


0.44 


0.56 


0.59 


Elastic+CEX 


1.088 


-0.142 


0.50 


0.65 


0.68 


All four channels 


1.103 


-0.203 


0.73 


0.94 


0.99 
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